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Description : 

L'invention concerne Hmagerie ultrasonore, par exemple, medicate. En particulier, 
Invention concerne un precede pour combiner des images pour former une image de 
combinaison en deux ou trois dimensions. Dans de tels precedes de combinaison spabale, 
5 I'objet image peut etre un plan, done en deux dimensions ou un volume, done en trois 
dimensions. 

Les brevets US 6,126,598 et US 6,390,980 decrivent des appareils et des precedes 
pour realiser des combinaisons spatiales dlmages ultrasonores en temps reel. La 
combinaison spatiale est une technique d'lmagerie selon laquelle les echos d'une able 
10 obtenus a partir de points multiples ou selon des directions dHnsonification multiples, sont 
melanges au sein d'une seule image dite de combinaison en combinant les donnees par 
exemple par moyennage, filtrage lineaire ou non-lineaire. Llmage de combinaison presente 
generalement un bruit (speckle en anglais) plus feible et une meilleure definition ou 
continuity des structures 2D ou 3D que les images ultrasonores conventionnelles produites a 
15 partir d'une unique direction d-insonification. En deux dimensions, plusieurs images sont . 

ainsi obtenues pour une meme position d'une barrette d'elements transducteurs en faisant j> 
varier I'angle du faisceau ultrasonore par rapport a la barrette. La combinaison damages est| 
alors equivalente a un moyennage temporel, les directions d'insonification etant sondees 
i'une apres 1'autre pour une meme position de*fcHSarrette, H 
20 De multiples techniques ont ete proposee/pour etendre le principe a des donnees en trois J 
dimensions. Parmi ces techniques, on peui citer la demande de brevet WO01/69282 qui * 
propose d'utiliser une acquisition dlmages ultrasonores par une barrette lineaire se ■ w 
deplacant par pas dans une direction perpendiculaire aux plans d'.mages. En multipliant les . ^ 
directions de deplacement de cette barrette ID, il devient alors possible d'acquerir plusieurs 
25 volumes ultrasonores imageant le meme milieu. Uobjet image est alors un volume. Apres 
avoir realigne les differents volumes entre eux, par exemple a partir des donnees d'un 
senseur de position, une combinaison de ces differents volumes est alors possible pour 
former un volume de combinaison. Compare a une seule acquisition d'un volume, ce volume 
• de combinaison offre un meilleur rapport signal sur bruit et augmente la definition et la 
30 continuite des structures 2D ou 3D du milieu observe. Bien entendu, a des fins de 
visualisation par exemple, on peut ne feire cette combinaison que dans un.plan bien 
particulier. La technique d'acquisition de donnees tri-dimensionnelles decrite ci-avant permet 
aisement dlmager le milieu dans un plan de coupe particulier. Par exemple une interpolation 
est realisee a partir de I'acquisition du volume echographique. Les differents volumes acquis 
35 par balayage du milieu avec differentes directions de deplacement de la barrette lineaire 
permettent done d'obtenir aisement differentes images du milieu dans un meme plan de 
coupe donne qui definit alors I'objet image et ainsi de construire une image de combinaison 



dans ce plan de coupe. Une operation simple et classique de combinaison spatiale consiste 
simplement a effectuer un moyennage de ces images. 

L'invention est liee aux considerations suivantes : 

Dans Tart anterieur WO01/69282, !es cibles dans un volume sont observees selon 
plusieurs directions, par exemple n directions. Des images de differentes resolutions sont 
alors combinees car I'image 2D ou 3D d'un refiecteur ponctuel (eest-a-dire la fonction de 
reponse du systeme d'imagerie a une impulsion) est fortement anisotrope. Cela provient des 
caracteristiques physiques de la sonde et du type dlnsonification, de focalisation et de 
formation de vote utiiisees. Du fait de la technique d'acquisition de donnees 3D decrite ci- 
avant, cette reponse du systeme est fortement dependante de la direction de deplacement 
de la barrette lineaire d'acquisition. Aussi un contour ou une partie de contour peut-il etre 
mieux visible ou untquement visible et blen delimite sur une image particuliere. 

Apres combinaison des differentes acquisitions du milieu avec plusieurs directions de 
deplacement de la sonde, I'image de combinaison presents une resolution plus isotrope ou 
encore une definition moyenne des contours car on synthetise ainsi une reponse moyenne 
du systeme. Cependant, la meilleure resolution ou la meilleure definition de contour 
presentee par au moins une des images combinees (par exemple, celle qui appartient a un 
des plans acquis par la barrette et qui a la ciuartiEe d'une image echographique 
conventionnelle) n'est pas conserves Cela est prejudiciable pour I'image de combinaison 
finalement obtenue. 

Un but de Invention est de profiter de la reduction du bruit en melangeant ies 
differentes Images tout en preservant la meilleure resolution possible et la plus grande 
definition ou continuite des structures obtenues possible. Ainsi, un but de (Invention est-il de 
permettre d'eviter de perdre Ies qualites d'une image particuliere au sein des images a 
combiner. 

En effet, un procede pour combiner des images conforme au paragraphe introductif 
est remarquable selon Pinvention en ce qui! inclut Ies etapes de : 

- recherche de contours representatifs d'une interface sur Ies images a combiner, ladite 
etape de recherche etant destinee a definir des zones d'interet au voisinage desdits contours 
representatifs, 

- analyse des zones d'interet, ladite etape d'analyse etant destinee a attribuer des poids aux 
points desdites zones d'interet et aux points correspondants auxdites zones d'interet sur Ies 
differentes images a combiner, 

- construction d'une image de combinaison, un point de I'image de combinaison 
correspondant a un point d'au moins une z<?ne d'interet 6tant obtenu a partir d'une 
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ponderation des points correspondants sur les images a combiner, ladite ponderation 
utilisant les poids attribues dans ladite etape d'analyse. 

Le resultat est que les meilleures donnees sont fiavorisees par les choix de ponderation en 
presence de structures anatomiques detectees comme contours. Ainsi, aucune structure 
5 anatomique n'est perdue ou trop attenuee. Hors des zones ou des structures sont presentes, 
Hmage est de bonne qualite car la combinaison classique permet de moyenner les images et 
done d'amenuiser le bruit qui est, lui aussi, moyenne. Le prindpe selon invention peut etre 
etendu a un nombre entier quelconque d'images. 

L'etape d'analyse inclut avantageusement des moyens de determination de caracteristiques 
10 des zones d'interet. Ces moyens de determination peuvent etre independants ou combines 
entre eux. Un certain nombre de modes de realisation de Hnvention correspondent a ces 
differents moyens de determination. 

Dans un premier mode de realisation, I'etape d'analyse utilise une detection des 
concordances des zones d'interet. Cette concordance est directement utilisee pour la 

15 determination des poids. Par exemple, cela permet qu'un contour detecte sur une unique , 
image et qui, par consequent, disparaitrait lors d'une moyenne des images, ne soit pas ^ 
affecte par la combinaison. Par exemple, les poids seront alors tous nuls pour les points 
correspondants sur les images distinctes de celle ou il a ete detecte. ^ \ 

Dans un mode de realisation avantageux, I'eJapetJ'analyse utilise une etude du contraste au 

20 sein des zones dlnteret En particulier les fcontrastes des contours sont ainsi evafues et 

utilfses pour calculer les poids. Les poids- pourront alors etre determines quantitativement en 
fonction de la valeur d'un facteur de confiance de la zone d'interet determine a partir de 
I'etude du contraste. 

Enfin, dans un mode de realisation prefer^, particulierement interessant pour une acquisition 
25 d'images r^alisee selon les principes presentes dans la demande de brevet WO01/69282, 
I'etape d'analyse utilise une etude de !a resolution. Pour une tefle acquisition, les images 
presentent en effet des resolutions differentes et il peut etre judideux de privilegier les 
images pour lesquelles la resolution est ta meilleure. 

Tout autre critere possible pour juger de la qualite d'un contour representatif d'une interface 
30 peut aussi etre choisi. On comprend bien id que la combinaison d'image prend le meilleur de 
chacune des images pour former Hmage globale de combinaison. En dehors des. structures 
anatomiques, le bruit est classiquement reduit par moyennage et, sur les structures, une 
meilleure impression visuelle est obtenue. 

L'invention peut done etre mise en ceuvre dans tout appareil destine a realiser une 
35 combinaison spatiale d'images d'un milieu. EUe concerne par consequent un dispositif 

destine a etre integre dans un appareil dlmagerie ultrasonore incluant des moyens pour 
realiser les etapes du precede precedemment presente. 
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Dans une de ses applications, 1'invention concerne done particulierement le domaine medical 
oil des images de grande precision sont les garantes d'un diagnostic juste et affine. 

L'invention sera mieux comprise a la lumiere de la description suivante de quelques 
5 modes de realisation, faite a titre d'exemple et en regard des dessins annexes, dans 
lesquels : 

la figure 1 est un diagramme schematique d'un appareii ultrasonore au sein duquel est 
avantageusement implementee l'invention, 

la figure 2 est un diagramme schematique d'un processeur destine a realiser un procede 

10 selon l'invention, 

la figure 3 est une illustration sur un exemple des effets d'un procede selon l'invention sur 

une combinaison dfmages, 

la figure 4 presente des exemples de techniques d'acquisition damages en deux dimensions 
permettant de reconstruire un volume, 
15 les figures 5a et 5b illustrent deux coupes orthogonales de.la fonction de reponse a un point 
ou fonction de reponse du systeme, une coupe etant faite dans le plan de la sonde et I'autre 
dans un plan orthogonal pris dans la direction de deplacement de la sonde, 
les figures 6a et 6b presentent deux types d'acquisition alternatives d'images a combiner 
selon un procede conforme a l'invention. 

La description qui va suivre est presentee pour permettre a un homme du metier de realiser 
et de faire usage de l'invention. Cette description est foumie dans le contexte de la demande 
de brevet et de ses exigences. Des alternatives diverses.au mode de realisation prefere 
seront evidences a lliomme du metier et les principes generiques de Invention exposes ici 
25 peuvent etre appliques a d'autres mises en oeuvre. 

Ainsi, la presente invention n'est pas censee etre limitee aux modes de realisation decrits 
mais plutot avoir la portee la plus large en accord avec les principes et les caracteristiques 
decrites ci-apres. 

Les remarques suivantes concernent les signes de reference. Des entites similaires sont 
30 designees par lettres identiques dans toutes les figures. Plusieurs entites similaires peuvent 
apparaitre dans une seule figure. Dans ce cas, un chiffre ou un suffixe est ajoute a la 
reference par lettres afin de distinguer des entites similaires. Le chiffre ou le suffixe peut 
etre omis pour des rafsons de convenance. Ceci s'applique pour la description ainsi que pour 
les revendications. 
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La figure 1 represente un systeme d'imagerie ultrasonore au sein duquel est 
avantageusement implementee la presente invention. Une sonde PRB incluant une barrette 
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ARR de transducteurs permet dMnsonifier un milieu. Une emission adequate et une formation 
de vote adaptee, comme decrit par exemple dans !e brevet WO01/69282 permettent de 
reconstruire un signal Radio Frequence (signal RF) pour differents angles d'insonification du 
milieu et d'obtenir ainsi des champs damages figures par un rectangle et des 
parallelogrammes en pointille. 

Trois groupes de lignes d'insonification sont indiques sur la figure et notes A, B & C, chaque 
groupe etant incline a un angle different par rapport a la sonde PRB. Remission des 
faisceaux est controlee par un emetteur TRM qui controle le decalage de phase et la duree 
d'acbvation de chaque element de la barrette ARR de transducteurs : ceci permet d'emettre 
chaque faisceau a une origine donnee sur la barrette de transducteurs et selon un angle 
predetermine. Les signaux regus par les elements de la barrette ARR de transducteurs sont 
digitalises par un convertisseur analogique/numerique et injectes dans un module digital de 
formation de vote DBF. Le module digital de formation de voie DBF retarde et somme les 
echos regus sur chaque element de la barrette pour former une sequence d'echantillons 
echographiques coherents et focalises le long de chaque ligne echographique. L'emetteur 
TRM et le module de formation de voie DBF sont pilotes par le controleur de systeme SYGy. 
repondant aux reglages fournis par un operateur via une interface utilisateur USI. Le module 
de controle du syst&me SYC controle l'emetteur TRM de maniere h transmettre le nombre^ 
desire de groupes de lignes d'insonification jjux^frgles souhaites dans les gammes d'enerfcjie 
et aux frequences desirees. Le controleur cJe systeme SYC controle aussi le module digital de 
formation de voie DBF de maniere a correctement retarder et combiner les echo requs pcmr 
les ouvertures et les profondeurs d'images utillsees. Les signaux d'echo sont flltres par unk 
filtre digital programmable FIL, qui definit une bande de frequences tfinteret puis transmis. 
au detecteur DET qui extrait I'enveloppe du signal et fournit I'arnplitude de cette derntere. 
Selon le principe de la presente invention, les signaux d'echos digitaux sont traites par 
combinaison spatiale dans un processeur P. Les signaux d'echo digitaux sont d'abord 
pretraites par un pre-processeur PRP. Le pni-processeur PRP peut eventuellement ponderer 
les echantillons de signaux avec un facteur de poids qui est une fopction du nombre 
d'images utilisees. Le pre-processeur peut aussi ponderer les points situes sur les bords des 
images superposees de maniere a lisser les transitions au niveau desquelles change le 
nombre d'echantillons ou d'images qui sont combinees. Les echantillons de signaux 
pretraites peuvent ensuifce subir un re-echantillonnage dans un re-echantillonneur RES. Le 
re-echantillonneur RES peut realigner dans Tespace les donnees tfune image ou les pixels 
d'une image sur I'espace d'affichage. Ensuite, les images sont combinees dans un module de 
combinaison CMB. La combinaison peut comprendre une somme, une moyenne, une 
detection de pic ou d'autres moyens de combinaison. L'invention est avantageusement mise 
en oeuvre au sein de ce module de combinaison CMB. Enfin un post-traitement est realise 
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par un post-processeur POP. Le post-processeur POP corrige la dynamique des valeurs 
combinees dans un intervalle de valeurs compatibles avec Paffichage. Le post-traitement 
peut etre racilement mis en beuvre a I'aide de look-up tables (LUT) et peut realiser 
simultanement la compression et la normalisation de maniere a ce que les valeurs 
5 combinees soient convenables pour 1'affichage. Pour assurer la compatibilite entre les tallies 
d'images issues du module POP et les tallies admissibles par le module d'affichage DIS et la 
carte video VIP, les images combinees sont re-interpolees par un scan-convertisseur SCC et, 
par exemple, stockees dans une memoire MEM. Apres scan-conversion, les images 
combinees de maniere spatiale sont traitees pour affichage par un processeur video VTP et 
10 affichees sur des moyens d'affichage d'images DIS. 

La figure 2 presente un mode de realisation d'un module de combinaison CMB permettant de 
realiser un precede selon Invention. Ce module de combinaison definit un dispositif destine 
a etre mis en ceuvre au sein d'un appareil dlmagerie ultrasonore selon Invention. Le 
module de combinaison CMB recoit les donnees de plusieursomages. Pour des raisons de 
simplification de la description, les donnees de deux images IMP], IM[J] seulement sont 
presentees en entree du module de combinaison CMB. Ces donnees ont avantageusement 
subi un premier pre-traitement classique tel que decrit ci-avant. Sur cette figure et les 
suivantes, les differentes notations designent generiquement, pour chaque image et pour 
chaque pixel ou voxel : les zones d'interet IA^frespondant a la presence ou non de 
structures dans les images, les poids W affectes a chaque pixel ou voxel de chaque image 
avant I'etape de combinaison generant fimage IM C . Les indices [i] et [J] permettent de 
distinguer les deux images entree du processus de combinaison CMB. 
La description de ce mode de realisation sera realisee pour un objet image de nature plane 
et done des images en deux dimensions. La figure 3 donne une illustration de deux images 
avantageusement traitees selon un procede de combinaison d'images selon (-invention. Les 
deux images representees sur la figure 3 peuvent resulter de differentes sortes d'acquisition 
d'une image du milieu. 

Dans une premiere application de Invention, les deux images IM[I] et IM[J] sont acquises 
pour une meme position de la barrette et pour differents angles d'insonification. La barrette 
ARR est alors positionnee selon le cas 1 represent sur la figure 3 par une barrette situee 
dans le plan de chacune des deux images. 

Dans une seconde application de ('invention, les deux images IMP] et IM[J] sont extraites de 
deux acquisitions tri-dimensionnelles comme decrit dans la demande de brevet 
WO01/69282 : elles correspondent alors a une meme coupe des deux volumes 
echographiques. En regie generate, ces images resultent d'une interpolation des donnees 
issues de chaque volume 3D et ne correspondent pas a une coupe acquise lors du 
deplacement de la barrette lineaire. Les deux Images (ou plus si on dispose de balayages 
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selon d'autres directions) correspondent au meme plan physique et peuvent done etre 
combinees selon les precedes decrits par exemple dans la demande de brevet WO01/69282. 
La figure 4 illustre le choix de tels plans de coupe dans des cas particuliers ou le plan de 
coupe est parallele ou orthogonal au plan de la sonde. Sur la figure 4, la barrette ARR etant 
orientee suivant la direction Y, la fleche AR designe la direction de balayage. Llmage BSC 
est ainsi acquise dans le plan de la barrette (Bscan en anglais). Llmage ESC est une coupe 
du volume balaye par le faisceau emis par la barrette ARR dans une direction orthogonale au 
plan de la barrette ARR (Elevation scan en anglais). Llmage CSC est une section orthogonale 
aux deux plans des images precedentes BSC et ESC indexee par la profondeur. 
Sur la figure 3, par souci de clarte et sans perte de generality IMP] etant par exemple une 
image BSC pour un angle determine dlnsonification, IMp] est alors par exemple une image 
ESC resultant d'une coupe pour un balayage dans une direction orthogonale a celle 
permettant I'acquisition de I'image IMffJ. La barrette ARR est alors posftionnee (et deplacee 
pour I'image IM[1]) selon le cas 2, e'est-a-dire, perpendiculaire au plan de I'image. 
Selon I'invention, un module de recherche de contour CNT recoit.les donnees damage. Le 
fonctionnement de ce module de recherche de contour QfT est avantageusement base sur 
les principes connus de recherche de contour correspondant a des transitions de niveaux du ^ 
signal echographique correspondant a des Interfaces significatives du milieu image. ^ 
L'homme du metier peut choisir une rechercbe»manuelle soft en detourant des zones 
d'interet, soft en dessinant des contours superposes sur des zones d'interet des images * 
echographiques. L'homme du metier peut'egalement choisir une recherche automatique on J, 
semi-automatique des contours. II peut s'atder de techniques d'aide au trace de contours, >. 
par exemple decrites de maniere generique dans W.A. Barrett et E.N- Mortensen, i, 
« Interactive Segmentation with Intelligent Scissors », Graphical Models and Image 
Processing, 60, pp 349-384, 1998. L'homme du metier a aussi a sa disposition tous les 
moyens de pre-segmentation d'images, d'extraction de contours... Parmi les moyens 
d'extraction de contour, on peut citer, par exemple, les filtres de Canny, de Shen, de 
Deriche... II y a ici une grande dlversite de choix possibles. Pour la realisation de I'invention, 
il est seulement necessaire d'obtenir en sortie du module de recherche de contour CWT des 
donnees permettant d'isoler un certain nombre de points de I'image detectes comme 
appartenant a des transitions significatives dans les images et correspondant a des 
interfaces presentes dans le milieu. Ces points permettent alors de definir au moins une 
zone d'interet IA dans leur voisinage. 

En effet et optionnellement, les contours C[I] et CP] constitues par les points detectes 
comme appartenant a des transitions peuvent etre grossis. Cela a pour but de supprimer les 
erreurs d'estimation dues au bruit et aux effets de pavage, de prendre en compte 
I'etalement spatial du a la reponse non parfaite du systeme d'acquisition, de considerer que 
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les interfaces ne sont pas infiniment fines. Dans ce cas, les zones d'interet IA sont definies 
par ces contours grossis. 

En se referant a la figure 3, les contours Cp] et Cp] sont recherches et les zones d'interet 
IA[I] et IA[J] sont definies sur les deux images IMP] et IMp]. Les zones d'interet IAp] et 
5 IAp] sont, sur la figure 3, le resultat d'un grossissement des contours detectes C[I] et Cp]. 
Ce grossissement est volontairement exagere. 

Les zones d'interet IAp], IAp] sont ensuite analysees dans un module d'analyse ANA. 
Dans un premier mode de realisation de I'invention, i'etape d'analyse ANA s'attache a 
determiner la concordance des presences de contours sur les differentes Images. Le but de 

10 I'etape d'analyse ANA est alors de determiner la superposition des zones d'interet sur les 
differentes images, Sur la figure 3 et pour simplifier ('explication, les zones d'interet 
resultant de cette recherche IAp] et IAp] (detectees sur les images IMp] et IMp]) sont 
representees en surimpression sur les images IMp] et IMP] (pointilles). On recherche done 
les points communs aux zones d'interet detectees sur chaque image, Lorsque des zones 

15 d'interet se superposent, le contour contenu dans ces zones d'interet a ete detecte sur les 
deux images et les deux zones d'interet peuvent done etre combinees en accordant a ces 
points des poids egaux. Les points appartenant a une zone d'interet detectee sur une image 
et pas a une zone cf interet sur I'autre image se voit attribuer un poids maximal sur I'image 
ou la zone d'interet a ete detectee et minfmal-pdur I'autre image. La ponderation est alors 

20 definie pixel par pixel de la zone d'interet. Ainsi, si on traite n images, 2 puissance n choix 

seront possibles pour chaque pixel. On peut par exemple decider d'affecter un poids presque 
nul en I'absence de zones d'interet et un poids de 1 dans le cas contraire. Par exemple, le 
tableau suivant pr^sente ('attribution des poids pour une combjnaison sur 2 images. Les 
deux premieres colonnes definissent par 0 et 1 1'appartenance d'un pixel a chacune des deux 

25 zones d'interet. W represente les poids et P les poids normalises, a etant un reel compris 
entre 0,5 et 1.. 



IA[I] 


IA[J 
] 


W[I] 


wpj 


W[I] + W[J] 


L J W[I] + W[J] 


1 


1 


0.5 < a <, 1 


0.5 < a <. 1 


0.5 


0.5 


1 


0 


0.5 < a <, 1 


0<l-a <0.5 


0.5 < en <1 
1 si a = 1 


0<l-a<0.5 
0 si a = 1 


0 


1 


0<l-a <0.5 


0.5 < a < \ 


0<l-a<0.5 
0 si a = l 


0.5<a<l 
1 si a = 1 


0 


0 


0$l-a <0.5 


0<l-a<0.5 


0.5 


0.5 



La valeur du pixel correspondant de I'image de combinaison est alors egale a : 
W\I] JM[/] + W[J)JM[J] = miM[I] + nJ] j M [J) 

c w[i]+w[j) 

Cette equation est generalisable a tout nombre n entier damages. La derniere ligne du 
tableau decrit en fait le cas classique des zones hors zones d'interet et done la combinaison 
5 spatiale damages classique. Une carte de poids peut ainsi etre obtenue pour les zones 
d'interet IA et pour les zones hors zones dlnteret notees IM-IA. La carte des poids peut 
done etre definie sur la totalite de I'image. 

Dans un mode de realisation particulierement avantageux, I'etape d'analyse ANA 
determine un indice de confiance pour un contour sur chacune des images ou il a ete 
10 detecte. Cela permet d'etablir une echelle entre les images en fonction de la qualite du 

contour. Ce mode de realisation utilise les resultats de I'etude de concordance des contours. 
Des poids non constants sur les points de la zone d'interet (par exemple autres que 0 et 1) ^ 
peuvent alors etre attribues en fonction, par exemple, du contraste observe sur les images a 
combiner ou de tout autre mesure de confiance. Ces poids'sont alors dependants de la i 
15 position du pixel ou du voxel considere. Ce contraste est determine a I'aide de moyens ?? 
classiques par calcul de grandeurs (telles le gradient, la derivee, par exemple) sur HmageS 
puis comparaison de ces grandeurs sur les differents points des zones d'interet concordanfes 
d'une image a I'autre. Dans I'exemple de la-ffgure 3, les points du contour qui appartienndnt 
aux deux zones dlnteret sont plus contrastes (figure par un trait plus large) sur I'image 
20 IMP]. Cette difference de contraste peut etre detectee par un calcul de gradient par * 
exemple puis une comparaison entre les deux images. L'appreciation du contraste peut atissi 
etre laissee a un utilisateur qui peut alors choisir une ou plusieurs images qu'il juge avoir fe 
meilleur contraste. Ainsi les points de I'image IM[J] se voient attribuer des poids plus 
importants sur la zone commune que les points correspondant de I'image IMP]. Le meilleur 
25 contraste observe sur le contour peut, par exemple, provenir d'une meilleure reflexion des 
ondes ultrasonores due a un angle d'emission du faisceau ultrasonore mieux adapte a 
I'interface que represente ledit contour. Cette interface est par exemple un changement de 
tissu dans le milieu observe. Enfin, afin de preserver une meilleure impression visueHe lors 
de I'observation de I'image de combinaison IMc, les images de poids normalises Pfl] et P[J] 
30 peuvent etre lissees avant le calcul de IMc de fagon a realiser une combinaison plus continue 
visuellement. 

Dans le mode de realisation prefere de Invention, I'etape d'analyse ANA met aussi en 
oeuvre une evaluation de la resolution sur 1es zones dlnteret des images a combiner. Ce 
mode de realisation utilise les resultats de I'etude de concordance des contours. Sur les 
35 differentes images ou des zones dlnteret communes ont ete detectees, une evaluation de la 
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resolution peut etre utilisee pour le calcul des poids. Cela est interessant lorsque les images 
a combiner ne presentent pas la meme resolution. Cest par exemple lecas pour deux 
images issues de la meme coupe de deux volumes echographiques acquis avec deux 
directions orthogonales de deplacement de la sonde ARR comme illustre sur la figure 4 : une 
5 des images correspond alors a une image de type BSC et est bien resolue, I'autre image 
correspondant a une image de type ESC plus grossiere. Dans le cas general, les images a 
combiner sont acquises avec une barrette lineaire de transducteurs orientee successivement 
selon plusieurs directions differentes lors des differentes acquisitions des volumes 
echographiques d'un meme milieu. Dans ce cas, la fonction de transfert resultant d'un echo 
10 d'un faisceau ultrasonore sur un point (reponse du systeme a une impulsion nommee Point 
Spread Function en anglais) est fortement anisotrope. 

Les figures 5a et 5b represented deux coupes de la fonction de transfert du systeme 
dlmagerie, vues dans deux plans de coupe particuliers. Pour simplifier les explications, on 
considere deux volumes echographiques acquis avec des deplacements de la barrette 

15 lineaire orthogonaux. On choisit d'observer deux memes ccjupes de chacun de ces deux 
volumes. La premiere coupe est perpendiculaire au deplacement de la sonde utilisee pour 
acquerir le premier volume. La seconde coupe est parallele au deplacement de la sonde 
utilisee pour acquerir le second volume. En se referant a la figure 3, la premiere coupe est 
une coupe de type BSC acquise dans le plan de^fa^onde (voir la figure 5a) et la seconde 

20 coupe est une coupe de type ESC acquise 4 dans un plan orthogonal au plan de la sonde (voir 
la figure 5b). Les deux figures 5a et 5b illustrent, pour un meme plan de coupe du milieu 
(mage, la resolution associee des images pour les differents points des plans etudies. Ces 
fonctions, resultant de I'echo sur un point, dependent de la profbndeur. Une degradation 
minimale d'echo E2 se sftue dans le plan de focalisation ou se situe un point 2, au centre de 

25 la tache. Les echos El et E3 sur des points 1 et 3, situes au centre des taches El et E3, 

presentent une resolution moins bonne que celle de I'echo E2 d'un point situe dans le plan 
de focalisation. On remarque d'autre part que ces fonctions qui represented le meme plan 
du milieu, mais issues de deux acquisitions du milieu avec des directions differentes de 
deplacement de la sonde, sont differentes. Elles illustrent la degradation des images 

30 acquises et interpolees par un balayage d'un milieu etudie a I'aide d'une barrette lineaire de 
transducteurs. La resolution laterale est en effet meilleure pour la figure 5a car elle 
correspond a la reponse du systeme dlmagerie dans le plan de la barrette. Cette resolution 
atteint sa degradation maximale dans le plan orthogonal au plan de la barrette, comme 
illustr£ sur la figure 5b. La combinaison de ces images melange done differentes resolutions 

35 spatiales. II en resulte une image presentant une resolution plus fsotrope. Cependant, la 
resolution globale est moins bonne que celle de la figure 5a. Selon I'invention, il peut etre 
ainsi avantageux de priviiegier I'image 5a acquise dans le plan de la barrette ARR par 
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rapport a celle 5b acquise par coupe dans une direction perpendiculaire a la barrette ARR. 
Ainsi une meilleure resolution globale sera obtenue. Ainsi, selon llnverrtion qui consiste a 
attribuer a des zones de certaines images une plus grande importance dans la combinaison, 
des poids plus importants sont avantageusement attribues aux zones d'interet concordantes 
5 detectees au sein damages acquises dans le plan de la barrette par rapport aux zones 
d'interet correspondantes detectees sur les images acquises par interpolation apres un 
balayage selon une direction parallele au plan de Timage. 

Les figures 6a et 6b decrivent le processus ci-dessus dans un cadre plus general si la 
combinaison est realisee a parb'r d'acquisitions faites avec des directions quelconques de 

10 deplacement de la barrette et dans un plan de coupe quelconque. Selon le mode de 

realisation prefere de finvention, est attribue un poids d'autant plus fort que la resolution est 
meilleure. Cela revient a attribuer un poids d'autant plus grand aux points d'une acquisition 
que celle-ci se trouve plus proche d'une coupe de type BSC En se referant a la figure 6 et 
pour un volume donne, on peut, par exemple, affecter un poids egal a la valeur absolue du 

15 cosinus de I'angle fait entre le plan de coupe de Timage PC et une coupe de type BSC 
Selon I'invention, on bent a preserver !e contraste sur les zones d'interet presentant des 
transitions et a effectuer un moyennage dans les parties du milieu observe ne contenant pas 
de structures anatomiques. Dans les parties ou les zones d'interet IA se superposent (voir , 
1'exemple de la figure 3), le contour contenu v dans"ces zones d'interet a ete detecte sur les 

20 deux images et les deux zones d'interet p^uvent done etre combinees en accordant a ces 

points des poids traduisant de la proximity de la coupe observee a une coupe de type BSC , 
Par exemple, le tableau suivant presente rattribution des poids pour une combinaison de 
deux images issues de deux volumes acquis selon les configurations decrites sur les figures i 
6a et 6b (acquisitions A et B). PC represente le plan de coupe a imager. 8 A est Tangle entre 

25 PC et un plan de coupe de type BSC pour I'acquisition A et G B est Tangle entre PC et un plan 
de coupe de type BSC pour Tacquisition &. Les deux premieres colonnes du tableau 
definissent par 0 et 1 I'appartenance d'un pixel a chacune des deux zones d'interet. W 
represente les poids et P les poids normalises, a etant un reel compris entre 0,5 et 1- 



IA[I] 


IAp] 


wp] 


W[J] 


1 J W[I]+W[J] 




1 


1 


ajcos^l : 

0.5 < a <1 


a.|cos^ y | 

0.5 < a < 1 


|cos^ 
Icos^l + lcos^j 
au lieu de 0.5 


|cos^j 
|cos^| + |cos^ y | 
au lieu de 0.5 


1 


0 


0.5 < a <1 


0<l-« <0.5 


0.5 < a < 1 
1 si a = 1 


0<l-a <0.5 
0 si a = 1 
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0<l-a <0.5 


0.5 <<* <1 


0<l-a <0.5 
0 si a —I 


0.5 < a < 1 
1 si a = 1 
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0 


0. 


0<\-a <0.5 




0.5 


0.5 








0<l-a <0.5 







La valeur du pixel correspondant de Hmage de combinaison est donnee par Inequation : 

Ainsi, cette ponderation permet de privilegier les donnees les mieux resolues, de ne pas 
5 perdre dlnformation sur les transitions et a nouveati de moyenner les zones sans structures. 
Enfin, afin de preserver une meilleure impression vjsuelle lors de 1'observation de Hmage de 
combinaison IMc, les images de pdids W[I] et W[J] ou les images de poids normalises Pp] 
et P[3] peuvent etre lissees avant le calcul de IMc de fagon a realiser une combinaison plus 
continue visuellement. 

10 En conclusion, 1'etape d'analyse ANA inclut done avantageusement des moyens de 
determination de la qualite d'un contour selon plusieurs criteres qui peuvent etre 
^dependants ou combines entre eux. Ces criteres comptent done, par exemple, une 
detection de la concordance de zones d'interet, une etude du contraste des contours, une 
etude de la resolution et tout autre critere possible pour juger de la qualite d'un contour 

15 representatif d'une interface. On comprend bien ici que la combinaison dlmage prend le 

meilleur de chacune des images pour formervlTmage globale IMc illustree sur la figure 3- Les 
ponderations sont realisees de maniere a 9e que les meilleures donnees soient favorisees en 
presence de structures anatomiques detectees comme contours, qu'aucune structure 
anatomique ne soit perdue et que, hors des zones ou des structures sont presentes, llmage 

20 soit de bonne qualite. En effet, prendre la moyenne des images sur les zones hors des zones 
d'interet permet de moyenner le bruit et done de I'amenulser. Le principe selon rinvention 
peut etre etendu a un nombre entier quelconque d'images. 

Dans une mise en ceuvre particuliere de invention, la modification des poids est offerte a 
Tutilisateur qui peut choisir de donner une importance plus grande a une image particuliere 
25 pour une zone d'interet donnee, ladite image particuliere etant jugee avoir !a meilleure 

qualite. Notamment, I'inten/ention d f un expert peut etre particulierement utile pour juger du 
meilleur contraste. Le resultat visuel peut alors etre evalue et affine par un utilisateur qui 
jugera de la qualite finale de I'image de combinaison et modifiera les poids en fonction de sa 
perception. 

30 Les cartes de poids W[I], W[J] finatement obtenues sont ensuite utilises dans un module de 
combinaison CC qui realise la sommation ponderee des differentes zones d'interet et des 
points correspondents auxdites zones d'interet lorsqu'elles ne sont pas concordantes d'une 
image a une autre. Hors des zones d'interet des differentes images a combiner, une 
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combinaison classique telle que connue de I'homme du metier est realisee. L/image de 
combinaison IMc est alors finalement obtenue en sortie du processeur.- • 
Avantageusement, des moyens d'amelioration de fimage peuvent etre mis en oeuvre de 
maniere independante a Invention (notamment et par exemple les techniques de 
« deblurring » et d' « enhancement » en anglais). 

L'invention a ete decrite pour des combinaisons d'images a deux dimensions. H a ete 
vu quels types d'images en deux dimensions peuvent etre combines (avec des orientations 
angulaires de faisceaux differentes et avec des orientations de balayage de la barrette 
differentes pour un meme plan damage). II est alors possible de reconstruire un volume en . 
trois dimensions a partir des images combinees en deux dimensions si celles ci sont 
disponibles suivant au moins deux directions de I'espace. Cea revient simplement a traiter 
plusieurs plans de coupe de fagon reguliere de maniere a generer un espace 3D reconstruit 
II est aussi possible de reconstruire des plans de coupe « a la volee », par exemple, 
seulement la coupe demandee par le module de visualisatioQ pilote par I'operateur. Dans ce 
cas, il est alors avantageux d'introdufre la resolution parmi les criteres de calcul des poids 
associes aux differentes images intersectees par ce plan de coupe. II est aussi envisageable 
d'appliquer les principes de Kinvention a une combinaison spatiale de volumes. Auquel cas, 
I'etape de recherche de contours detecte des surfaces dans les volumes et I'etape d'analyse 
pondere des points de zones dlnteret qui sonfcdes-volumes dlnteret. Les principes selon 
Invention peuvent done etre generalises aTd'autres types de combinaison d'images et a 
d'autres types d'images. 

I! existe de nombreuses fegons de mettre en oeuvre les fonctions presentees dans les 
moyens et les etapes du procede selon l'invention par des moyens logiciels et/ou materiels 
accessibles a I'homme du metier. Cest pourquoi les figures sont schematiques. Ainsi, bien 
que les figures montrent differentes fonctions realisees par differents blocs, cela n'exclut pas 
qu'un seul moyen logiciel et/ou materiel permette de realiser plusieurs fonctions. Cela 
n'exclut pas non plus qu'une combinaison de moyens logiciels et/ou materiels permette de 
realiser une fonction. 

Bien que cette invention ait ete decrite en accord avec les modes de realisation 
presentes, un homme du metier reconnartra immediatement qu'il existe des variantes aux 
modes de realisation presentes et que ces variantes restent dans I'esprit et sous la portee de 
la presente invention. Ainsi, de nombreuses modifications peuvent etre realisees par un 
homme du metier sans pour autant s'exclure de I'esprit et de la portee definis par les 
revendications suivantes. 
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Revendications : 

1. Procede pour combiner des images d'un meme objet incluant les etapes de : 

- recherche de contours representatifs d'une interface sur les images a combiner, ladite 
5 etape etant destinee a definir des zones d'interet au voisinage desdits contours 

representatifs, 

- analyse des zones d'interet, ladite etape d'analyse etant destinee a attribuer des poids aux 
points desdites zones d'interet et aux points correspondents auxdites zones d'interet sur les 
differentes images, 

10 - construction d'une image de combinaison, un point de llmage de combinaison 

correspondent a un point d'au moins une zone d'interet etant obtenu a partir d'une 
ponderation des points correspondants sur les images a combiner selon les poids attribues 
dans ladite etape d'analyse. 

». 

15 2. Procede selon la revendfcation 1, caracterise en ce que I'etape d'analyse met en oeuvre 

une etape devaluation d'une concordance des zones d'interet sur les images a combiner, les 
poids etant attribues aux differents points desdites zones d'interet et a leurs points 
correspondants en fonction de ladite concordance. 

20 3. Procede selon la revendication 2, caracterise en ce que I'etape d'analyse met en oeuvre 
une etape d'estimation du contraste au sein d'au moins deux zones d'interet presentes et 
concordantes sur deux images, les poids etant attribues aux differents points desdites zones 
d'interet en fonction dudit contraste estime. 

25 4. Procede selon la revendication 2, caracterise en ce qu'au moins deux images a combiner 
presentent une resolution differente et en ce que I'etape d'analyse met en oeuvre une etape 
devaluation de ces resolutions au sein d'au moins deux zones d'interet presentes et 
concordantes sur deux dites images, les poids etant attribues aux differents points desdites 
zones dlnteret sur lesdites deux images en fonction desdites resolutions. 

30 

5. Dispositif destine a etre integre au sein d'un appareil d'imagerie ultrasonore et destine a 
former une image de combinaison, resultat d'une combinaison spatiale d'images d'un meme 
objet, incluant des moyens de : 

- recherche de contours representatifs d'une interface sur les images a combiner, ladite 
35 etape etant destinee a definir des zones dlnteret au voisinage desdits contours 

representatifs, 
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- analyse des zones d'interet, ladlte etape d'analyse etant destinee a attribuer des poids aux 
points desdites zones d'interet et aux points correspondants auxdites zones d'interet sur les 
differentes images, 

- construction d'une image de combinaison, un point de runage de combinaison 

correspondant a un point d'au moins une zone d'interet etant obtenu a partir d'une 
ponderation des points correspondants sur les images a combiner selon les poids 
attribues dans ladite etape d'analyse. 

6. Dispositif selon la revindication 5, caracterise en ce que les moyens d'analyse mettent en 
ceuvre des moyens devaluation d'une concordance des zones d'interet sur les images a 
combiner, les poids etant attribues aux differents points desdites zones d'interet et a leurs 
points correspondants en fonction de ladite concordance. 

7. Dispositif selon la revendication 6, caracterise en ce que les moyens d'analyse mettent en 
ceuvre des moyens d'estimation du contraste au sein d'au moins deux zones d'interet 
presentes et concordantes surdeux images, les poids etant attribues aux differents points 
desdites zones, tifhtlret en fonction dudit contraste esb'me. 

8. Dispositif selon la revehdfcation 6, caracte^T=fn ce qu'au moins deux images a combiner 
presentent une resolution differente et en .'ce que les moyens d'analyse mettent en ceuvre 
des moyens d'evaluation de ces resolutions au sein d'au moins deux zones dlnteret 
presentes et concordantes sur deux dites images, les poids etant attribues aux differents - 
points desdites zones dlnteret sur lesdites deux images en fonction desdites resolutions. 

9. Appareil dlmagerie ultrasonore incluant un dispositif destine a former une image de 
combinaison, resultat d'une combinaison spatiale d'images d'un meme objet, selon I'une des 
revendications 5 a 8. 

10. Produit programme d'ordinateur comprenant des instructions pour mettre en oeuvre un 
procede pour combiner des images d'un meme objet selon I'une des revendications 1 a 4. 
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